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ABSTRACT 

Aims. The accretion-powered pulsar 4U 1626-67 experienced a new torque reversal at the beginning of 2008, after about 18 years of 
steadily spinning down. The main goal of the present work is to study this recent torque reversal that occurred in 2008 February. 
Methods. We present a spectral analysis of this source using two pointed observations performed by Suzaku in 2006 March and in 
2010 September. 

Results. We confirm with Suzaku the presence of a strong emission-line complex centered on 1 kcV, with the strongest line being the 
hydrogen-like Ne Ly a at 1 .025(3) keV. We were able to resolve this complex with up to seven emission lines. A dramatic increase of 
the intensity of the Ne Lyar line after the 2008 torque reversal occurred, with the equivalent width of this line reaching almost the same 
value measured by ASCA in 1993. We also report on the detection of a cyclotron line feature centered at ~37 keV. In spite of the fact 
that an increase of the X-ray luminosity (0.5-100 keV) of a factor of ~2.8 occuned between these two observations, no significant 
change in the energy of the cyclotron line feature was observed. However, the intensity of the -1 keV line complex increased by an 
overall factor of ~8. 

Conclusions . Our results favor a scenario in which the neutron star in 4U 1626-67 accretes material from a geometrically thin disk 
during both the spin-up and spin-down phases. 
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1. Introduction 

4U 1626-67 is a Roche-lobe filling binary system, where mat- 
ter from the companion flows from the LI point into an accre- 
tion disk (Reynolds et al., 1997). This low mass X-ray binary 
(LMXB) consists of a 7.66 s X-ray pulsar accreting from an ex- 
tremely low mass companion ((0.04 Mo for i = 18°, Levine 
et al., 1988). Although orbital motion has never been detected 
in the X-ray data, pulsed optical emission reprocessed on the 
surface of the secondary revealed the 42 min orbital period 
(Middleditch et al., 1981), confirmed by Chakrabarty (1998). 
Most likely the binary system contains a hydrogen-depleted 
secondary to reach such a short orbital period (Paczynski & 
Sienkiewicz, 1981). The faint optical counterpart (KZ TrA, 
VM7.5) has a strong UV excess and high optical pulse fraction 
(McClintock et a|., 1977, 1980). 

The sign of the torque in 4U 1626-67 has reversed in only 
two occasions since its discovery by Uhuru (Giacconi et al., 
1972). Monitoring of the source by the Burst and Transient 
Source Experiment (BATSE) on board the Compton Gamma 
Ray Observatory ( CGRO ) starting in April 1991 , found the pul- 
sar spinning down, implying a changed sign in the accretion 
torque (Wilson et al., 1993; Bildsten et al., 1994). It was es- 
timated that the reversal occurred in mid-1990. More recently, 
our daily monitoring of 4U 1626-67 with Fermi/GBM starting 
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in 2008 August found the pulsar again spinning-up rather than 
spinning-down. Swift/BAiY observations allowed us to cover the 
evolution of this second torque reversal. The transition took 
place at around MJD 54500 (2008 Feb 04) and lasted approx- 
imately 150 days (Camero-Arranz et al., 2010). 

A 48 mHz quasi-periodic oscillation (QPO) has been de- 
tected in the X-ray emission (Shinoda et al., 1990; Kommers 
et al., 1998). Orland ini et al. (1998) inferred a neutron star mag- 
netic field in the range of (2 .4-6 .3 )x 10 12 G. To compute this 
magnetic field range a source distance of 5-13 kpc was assumed 
(Chakrabarty, 1998, and references therein). A -37 keV absorp- 
tion cyclotron feature was found in the 0.1-200 keV BeppoSAX 
spectrum (Orlandini et al., 1998). A more detailed study of this 
feature with Suzaku can be found in Iwakiri et al. (20 12). The X- 
ray broad-band continuum was fitted with low-energy absorp- 
tion, a blackbody, a power law and a high energy cutoff. 

Angelini et al. (1995) using the ASCA observatory discov- 
ered for the first time a strong emission-line complex centered 
on 1 keV in the X-ray spectrum of 4U 1626-67. The strongest 
line was identified as hydrogen-like Ne Lycr at 1.008keV. The 
strength of the neon emission compared to the expected iron L 
complex implied a large neon overabundance. 

We present a spectral analysis of 4U 1626-67 with Suzaku, 
using two observations from 2006 March and 2010 September. 
This allows us to better characterize the recent torque reversal 
that occurred in 2008 February. 
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Table 1. Best Broad Band Fit Spectral Parameters* 


Observation 

2V H 

a? 

norm 


E d 

raid 

^CRF 

Tfifloi* 

nOIIn lfl 

Flux/ 

constant 

factor'’ 

*?(DOF) 

400015010 
(Mar 2006) 

146(2) 

0.846(9) 

0.010(1) 

18.9(5) 

13(1) 

36.9^| 

3.6(6) 

0 229 +001 

-0.009 

127± 13 

0.55(9) 

0.85(2) 

0.8(fixed) 

138(205) 

405044010 
(Sep 2010) 

0.6(2) 

0.981(3) 

0.026(1) 

21.6(7) 

15(2) 

37.8(8) 

4.3(6) 

0.511(8) 

412!!^ 

1.58(4) 

1.21(4) 

0.8(5) 

1.20(206) 


fl (PHABS (8 GAUSSI ANS +GAUSSI AN+ BBODYRAD) +P0W)*H1GHECUT *GABS)*CON$T 
b x 10 21 cm -2 

c Photon Index with normalization units in photons keV’cm -2 ® -1 at 1 keV. 
d keV 

** The values for the cyclotron line depths are I4±4 and 18±4 for the Mar 2006 and Sep 2010 observations, respectively. 

' with normal =R| flfon /d^ 0 , where dio is the distance in units of 10 kpc, varying from 0.5-1.3 (assuming the 5 -13 kpc distance range reported for this source). 

I xlO -9 erg cm -2 s^ 1 (0.5-100 keV), based on the XIS 03 flux calibration; XIS 03 is a combined spectrum from XIS 0 and XIS3. 

3 Intercalibration factors (P1N/XIS 03 and GSO/XIS 03). The GSO/X1S03 factor was fixed to the » alue obtained for the second obsen alion since it could not be well 
constrained. 


2. Observations and Analysis 

The instruments on Suzaku (Mitsuda et al., 2007) are the X-ray 
Imaging Spectrometer (XIS; Koyama et al., 2007) CCD detec- 
tor covering the 0.3-10 keV band, and the Hard X-ray Detector 
(HXD; Takahashi et al., 2007) comprised of the PIN diode de- 
tector (PIN) covering the 10-70 keV band and the gadolinium 
silicate crystal detector (GSO) covering the 60-600 keV band. 
The XIS has four separate detectors, XIS 0-3, with XIS 1 being 
a backside illuminated CCD and the others frontside. XIS 2 was 
lost due to a micrometeor hit in late 2006. 

The first Suzaku 4U 1626-67 observation occurred in 2006 
March (MJD 53797) and thus before the torque reversal in 2008. 
During this observation the spacecraft aimpoint was optimized 
for the XIS detectors, and Suzaku was run in a 1/8 window mode 
with 1 s integration time, an effective area of 128 X 1024 pixels 
and an effective exposure time of -72.2 ks. The second observa- 
tion in 2010 September (MJD 55444) was also centered on the 
XIS aimpoint but run in a 1/4 window mode, with 2 s integration 
time, an effective area of 256 x 1024 pixels and an effective ex- 
posure time of -14.01 ks. For this observation only 3 XIS detec- 
tors were available, as we discussed in the previous paragraph. 
The Suzaku data were reduced with tools from the HEASOFT 
v6.11 package (Amaud, 1996) and CALDB release 20110804. 
Following the Suzaku ABC guide 1 , the first step of our analysis 
was to reprocess the data by running the FTOOL aepipeline for 
the XIS, the PIN and the GSO. The pipeline performs calibration 
as well as data screening. 

Due to thermal distorsion of the Suzaku spacecraft the opti- 
cal bench shows wobbling. The source’s PSF image is not cen- 
tered on a fixed position on the XIS CCDs but appears blurred. 
As part of the standard data reduction this effect is partially cor- 
rected by adjusting the spacecraft attitude file. The attitude file 
can be further improved by applying the FTOOL aeattcor2: 
based on the algorithm described by Nowak et al. (2011) the 
mean event positions as a function of time are calculated and 
used to further adjust the attitude. While we obtain sharp images 
for our first Suzaku observation, the second one suffers from 
additional attitude variability that can be present in Suzaku ob- 
servations since a change in the attitude control system on 2009 
Dec 18. In this case the FTOOL cannot fully correct the character- 
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istic "double images". According to the documentation 2 , how- 
ever, and since conservative pileup exclusion regions were ap- 
plied (see below), the XIS spectra are not significantly affected. 

We used the FTOOL pileest to estimate the amount of 
pileup in the XIS images and disregarded circular regions with 
> 8% pileup fraction during the spectral extraction. This is a 
careful approach since due to the asymmetric shape of the PSF 
only few pixels outside of the exclusion region reach pileup frac- 
tions above 2 - 4%. For the first observation the pileup fraction 
remained below this limit in the entire source extraction region (a 
circle of - 70" radius centered on the source image). In the sec- 
ond observation pileup fractions of up to 15% were reached in 
the center of the source image so the extraction region consisted 
of a circle of - 130" radius minus a central circular exclusion 
region with a radius of - 35" . 

The primary tool for extracting data products (spectra, 
lightcurves) from XIS data is xselect , which is part of the gen- 
eral HEAsoft distribution. Events in the XIS detectors were in 
the 3x3 and 5x5 editing modes. With xselect we applied to 
the cleaned events filters with the good times intenals avoid- 
ing telemetry saturation. Source and background regions were 
created from the image. The filtered events were then used to 
create source and background spectra for the desired extraction 
regions. Products from the two editing modes were added. The 
XIS spectra were grouped close to the energy resolution of the 
detectors (see Nowak et al ., 201 1). We created response matrices 
and effective area files with the xisrmfgen and jc issimarfgen tools, 
respectively. The calibration of the XIS is uncertain around the 
Si K edge, thus events in the 1. 8-2.0 keV range were excluded 
in the fitting. For more information about this issue we refer the 
reader to the "XIS Calibration around the Si K edge" web page 3 
(see also Iwakiri et al., 2012; Nowak et al., 2011; Suchy et al., 
2011, for different treatments on this matter). 

For the PIN, we extracted spectra from the cleaned event files 
using hxdpinxbpi. This tool produces the dead time corrected 
PIN source spectrum as well as the PIN background spectrum. 

The non X-ray background events provided by the Suzaku team 
can be downloaded from the HEASARC archive 4 . The PIN spec- 
tra were grouped to have a signal-to-noise ratio £ 10 in each 
energy bin, and we considered spectra between 14-70 keV. We 

2 ftp://legacy.gsfcnasa.gov/suzaku/doc/gencral/suzakumemo-2010- 
05.pdf 

3 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/sicaI html 

4 ftp://legacy.gsfc nasa.gov/suzaku/data/background/pinnxb_ ver2.0 .tuned 
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Fig.l. The main panels show the Suzaku broad-band spectrum before (left panel) and after (right panel) the 2008 torque reversal. Blue points 
denote XIS data (0.5-10 keV), PIN data points are in red (14-70 kcV) and GSO in purple (50-70 keV and 50- 1 00 keV, respectively). The residuals 
after fitting only the continuum model are shown in the middle panel, with no line emission and cyclotron line feature included. The bottom panel 
shows the residuals after fitting the line complex at -1 kcV, the -6.5 keV iron line and the -37 keV cyclotron line (see Tabic 1 for a list of fitted 
parameters). A zoom of the region of the line complex at -1 keV is also shown for clarity. 


point out that due to the changes in instrumental settings the PIN 
response matrices, from the CALDB database, are different de- 
pending on the epoch of the observation. 

The GSO spectra were created starting with the clean e\ ent 
files using the hxdgsoxbpi tool. The background events were 
downloaded from HEASARC 5 . Response files were taken from 
the CALDB database and the GSO correction ARF file was 
provided by the HXD instrument team 6 . The distributed GSO 
Non X-ray Background is grouped, so GSO spectra have to be 
grouped accordingly, using the spectral binning given by the file 
gsobgd64bins.dat available at HEASARC 7 . The GSO spectrum 
was restricted to the 50-70 keV range for the first observation, 
and to the 50-100keV for the second observation. 


3. Results 

The -0.5-100 keV broad-band spectrum obtained using XIS (a 
combination of detectors 0 and 3), PIN and GSO were fitted in 
XSPEC 12.7.0. We used a canonical model which allowed us 
to compare our spectral study with previous works (Angelini 
et al., 1995; Pravdo et al., 1979; Orlandini et al., 1998; Krauss 
et al., 2007; Jain et al., 2010, among others) and update the long- 
term X-ray flux history of 4U 1626-67 relative to the flux mea- 

5 ftp://legacy.gsfc .nasa.gov/suzaku/data/background/gsonxb_ vcr2.5 

6 http://www.astro.isas .ac.jp/suzaku/analysis/hxd/gsoarf/ 

7 ftp://legacy.gsfcjiasa.gov/docs/suzaku/analysis/gsobgd64bins.dat 


sured by HEAO 1 (Chakrabarty et al., 1997; Krauss et al., 2007; 
Camero-Arranz et al., 2010). 

The continuum model fitted to both observations in- 
cludes low-energy absorption, a blackbody component, and a 
power law with a high-energy cutoff at -20 keV (PHABS 
(BBODYRAD+POWLAW) HIGHECUT)*CONST. The con- 
stant factor (CONST) is used to account for intercalibration be- 
tween XIS, PIN, and GSO. The abundance was set to wilm 
(Wilms et al. 2000). Spectral fits for both observations are shown 
in Figure 1 . A broad iron line near 64 keV (GAUSS) and a 
-37 keV cyclotron feature (GABS) improved the fit for both ob- 
servations (see also Iwakiri et al ., 201 2). Table 1 presents the fit- 
ted values for these broad-band parameters from the final model 
(which included the low-energy lines described below). 

Initial fitting results for the broad band model for both obser- 
vations indicated evidence for a strong line complex at ~1 keV, 
shown in Figure 1, middle panels. Following Angelini et al. 
(1995), we added seven narrow Gaussian lines (listed in Table 2) 
to our model. The line widths were fixed at the ASCA values 
(JO eV for all lines except the 1362 keV line width was fixed 
at 57 eV). Table 2 lists the final fitted line parameters for both 
observations. Only five of the seven lines detected with ASCA 
were detected in the 2006 Suzaku observation, while six ASCA 
lines were detected in the 2010 Suzaku observation. For the tree 
ASCA lines not well detected in 2006 and the one in 2010 the 
intensities were fixed, since the errors could not be well con- 
strained. 
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Fig. 2. The top panel shows the 8 Gaussian lines from the -1 keV complex detected by Suzaku , before (left panel) and after (right panel) the 2008 
torque reversal. Here the dotted lines in the model denote the components that were fixed. The residuals after fitting only the continuum model are 
shown in the middle panel, with no line emission included. The bottom panel shows the residuals after fitting the line complex with those 8 lines 
(see Table 2 for more information about these lines). 


Table 2. Fit Results on Line Emission 0 


Observation 

Expected 

(keV) 

line Energy (keV) 

Mar 2006 Sep 2010 

A* 2 

Mar 2006 Sep 2010 
(dof =205) (dof=206) 

Intensity* 
Mar 2006 

Sep 2010 

EW(eV) 

Mar 2006 Sep 2010 

OHea 

0.568 

0.568(fixed) 

0.568(fixed) 

- 

- 

0.09(fixed) 

0.09(fixed) 

3.2(2) 

13(2) 

O Lyor 

0.653 

0.661(8) 

0.66(2) 

59 

5.4" 

3-0 

4(4) 

13(4) 

7(3) 

Ne Hear'* 

0.914 

0.930(8) 

0.929(8) 

105 

119 

1-3(3) 

7(1) 

8(1) 

14(1) 

Ne Lyar 

1.021 

1.025(3) 

1.031(2) 

878 

1699 

33(2) 

24(1) 

22.6(2) 

50.9(8) 

NeHc^ 

1.084 

1.084{fixed) 

1.13(3) 

— 

39 

0.2(fixed) 

3(1) 

3.7(2) 

6(2) 

NcLyfi 

1.210 

1.210(fixed) 

1.25(1) 

6 

96 

02(fixed) 

2(1) 

1.8(2) 

9(4) 

Mg Hear 

134 

132(2) 

138(2) 

36 

31 

0.4(1) 

1.8*7 

4(1) 

5(1) 

Ne Recombo* 

1362 

1.48(2) 

132(2) 

137 

154 

0.9(2) 

5(1) 

11.1(3) 

15(5) 

Fe K-fiuor. lincf 6.4 

63(2) 

6.78(7) 

59 

30 

0.8^ 

13(5) 

422(3) 

33(1) 


* the widths of the low -energy lines were fixed to the ASCA values (10 eV for all but the 1.362 keV line with 57 eV) 

* photons lO^cnT 2 *” 1 

c we note that, for the 2010 observation, the Ax 2 for this line is not well computed due to the spectral model not being able to adequately fit this line (sec Fig. 2, right panel). 
d the abundances of the lines at 0.914 and 1.084 keV can also be consistent with a blend of several L-shell iron line emission (see Angelini et al. 1995). 
r the Mg Ha line emission occurs at 1.471 keV (here not resolved, and thus not detected, by Suzaku), and therefore this Mg line may be contributing to this feature. 

■f for the iron line the widths were fitted (0.5(3) keV and 0.14(7) keV, respectively) 


The broad feature at 1 .362 keV could alternatively be ex- 
plained by a recombination edge. Fits to both observations using 
a model that replaced the 1 .362 keV Gaussian line with a re- 
combination edge yielded similar tfy 2 to our model containing 
the 1 .362 keV line. The edge (and temperature) obtained for the 
2006 and 2010 observations are 1.45±0.05 keV (42 eV) and 
1.47±0.03 keV (45^ eV) (fitting with a Gaussian line gives 
1 .48(2) and 1 .51(2) keV, see Table 2), respectively. A posteriori 
inspection of the residuals for both observations, suggested the 
presence of an extra Gaussian component at energy ~1.34keV. 
Therefore, up to 8 Gaussian components were included in the 
fit (see Fig. 2; A f of 283 for 205 dof, and 248 for 206 dof, 
for the first and second observations, respectively). Final fit pa- 


rameters for this line at ~1.34keV, tentatively identified as the 
He-like Mg, are listed in Table 2. The Ax 2 obtained by setting 
the flux of a particular line to zero in the best-fit model indi- 
cates the significance of each line. In addition, we fit the broad- 
band spectrum using different combinations of the XIS detectors 
yielding similar results. Other models accounting for emission 
from a hot diffuse gas, with line emissions coming from several 
elements, were also fit (meka and vmeka in XSPEC notation). 
We changed the elemental abundances while performing the fits 
with the meka and vmeka models, however they did not give rea- 
sonable fits. 

A luminosity increase occurred between the two observa- 
tions, changing from L<o.5-io*«v) «2.4xl0 34 d^ r to 6xl0 34 d^ 
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erg s" 1 . However, the value of the center of energy of the cy- 
clotron line did not exhibit a significant change. This behavior 
resembles the case of the Be/X-ray binary system A 0535+26 
• Caballero et al., 2009), but is contrary to sources like HerX-1 
and V 0332+53 (Pottschmidt et al ., 201 1 , and references therein) 
or GX 304-1 (Klochkov et al., 2012), which show positive or 
negative L-Ecrf correlations (see also Becker et al., 2012). 
On the other hand, the column density, A H , dropped by a fac- 
tor of ~3 while the black body temperature increased by a fac- 
tor of ~2. The rest of the broadband parameters were consis- 
tent between the two observations. Absorbed fluxes in the 0.5- 
'.OkeV, 2-10keV, 2-20 keV and 2-50keV bands are 0.18(4), 
0.159(2), 0.295(3), 0.51(3) xl0“ 9 erg cnrV 1 and 0.49(3), 
0.38(9), 0.64(1), 1.3(9) xlO" 9 erg cnrV 1 , for the 2006 and 
2010 Suzaku observations, respectively. 


4. Discussion and Conclusions 

4. t. Spectral Continuum 

Our results show that two years after the torque reversal in 2008 
the luminosity has reached almost the same level as in 1977 (see 
Fig. 3), the black body temperature has increased to ~ 0.6 keV 
and the pulsar shows the same spin-up rate. The system’s appear- 
ance in this state can be explained in terms of accretion onto the 
neutron star from a geometrically thin Keplerian disk provided 
the mass-transfer rate is limited to Af > 10 16 gs' 1 (Chakrabarty 
et al., 1997). 

The pulsar’s transition from spin-up to spin -down in 1990 
•vas accompanied by significant changes of spectral parame- 
ters. The black body temperature decreased from ~ 0.6 keV to 

- 0.3 keV and the power law photon index changed from ~1.5 
< Vaughan & Kitamoto, 1997) to ~0.6 (Vaughan & Kitamoto, 
1997; Angelini et al., 1995; Owens et al., 1997). The transition 
period lasted about 150 days and the source luminosity during 
this time was only a factor of 2 smaller than its highest level 
observed in 1977 (see Fig. 3). This indicates that the mass- 
transfer rate during the torque reversal in 1990 was in excess 
cf 5 x 10 15 g s _1 . 

Our Suzaku observations show that the pulsar’s behavior 
during the torque reversal in 2008 (from spin-down to spin- 
up phase) was different. While the transition period also lasted 
about 150 days and the black body temperature changed from 

- 0.23 keV back to ~ 0.6 keV in 2010, the photon index re- 
mained almost unchanged. The value of the order of unity ob- 
served in 2006 and 2010 is the same as the photon index ob- 
served by several missions before the torque reversal during the 
spin-down phase, e.g. ASCA (Angelini et al., 1995), BeppoSAX 
(Orlandini et al., 1998), Chandra and XMM (Krauss et al., 2007; 
Schulz et al., 2001). The torque reversal in 2008 occurred al 
slightly lower luminosity level than in 1990 and was followed 
by a brightening of the X-ray source on a time scale of more 
than 2 years. 

As we have shown in a previous paper (Camero-Arranz et al., 
2010), the torque reversals in 4U 1626-67 are unlikely to be as- 
sociated with alternation of the accretion flow geometry between 
the Keplerian geometrically thin disk (during the spin-up phase) 
and a spherical flow or a hot, geometrically thick disk (during the 
spin-down phase). The corresponding scenarios (see, e.g. Yi & 
Vishniac, 1999) encounter major difficulties explaining 150day 
transition period (see also Wijers & Pringle, 1999) and the low 
level of X-ray flux measured before and after the torque reversal 
in 2008. Furthermore, the mass-transfer rate during both of the 
observed torque reversals (> 5xl0 15 g s“ l ) was above the critical 


mass-transfer rate at which the transition of the accretion disk to 
the hot, geometrically thick state can be expected (Narayan & Yi, 
1995). Finally, observations of the strong emission-line complex 
centered on 1 keV and the iron K-fluorescence line at 6.4 keV 
during both spin-up and spin-down phases indicate that the tem- 
perature of the material surrounding the neutron star does not 
exceed a few keV. Conversely, the temperature of the material 
in a thick disk is close to the adiabatic temperature, which for 
the parameters of 4U 1626-67 is in excess of 10 keV. Thus, our 
results favor a scenario in which the neutron star in 4U 1626-67 
accretes material from a geometrically thin disk during both the 
spin-up and spin-down phases. 

The magnetospheric radius of a neutron star accreting ma- 
terial from a geometrically thin disk can be evaluated as r m ^ 

*(/i 2 /M(2GM ns ) 1/2 ) 1 , where p = (l/2)fU?^ s is the dipole 
magnetic moment, B> the surface field and the radius of 
the neutron star, k is a parameter accounting for the geome- 
try of the accretion flow, which in the case of disk accretion 
is limited to k 0.5 (Ghosh & Lamb, 1979; Wang, 1996; 
Bozzo et al., 2009). For the star to be in the accretor state its 
magnetospheric radius should not exceed the corotation radius, 

?cof = (GA/nsFf/4^ 2 ) 13 . Otherwise, the centrifugal barrier at the 
magnetospheric boundary would prevent the accretion flow from 
reaching the stellar surface (Shvartsman, 1970). This condition 
in the case of 4U 1626-67 is satisfied if the mass-transfer rate is 
M > M m io> where 

Mmin - 6 X 10 14 g S' 1 4 / >" S/3 ^.7 % 0) 

Here kq . 5 = k/ 0.5 is normalized to its minimum value, m is the 
mass and P7.7 the spin period of the neutron star in units of 
1.4 Mq and 7.7 s. The dipole magnetic moment of the neutron 
star, /i3o, is normalized to l0 3O Gcm 3 . According to our mea- 
surements of the energy of the cyclotron line feature with Suzaku 
(36.8 keV) the surface field of the neutron star in 4U 1626- 
67 is r 3 x 10 l2 G. This is consistent with previous estimate 
of the field strength made by Orlandini et al. (1998) from ob- 
servations of 4U 1626-67 with BeppoSAX. Putting this value to 
Eq. (1) and taking into account that 4U 1626-67 is a persistent 
pulsar one finds that the mass-transfer rate has never decreased 
below 10 15 g s~ l in the time period between 1977 and 2010. The 
rate of mass- transfer, therefore, may vary in the system up to 
an order of magnitude. It is remarkable that the observed X- 
ray flux (and, correspondingly, the mass accretion rate onto 
the surface of the neutron star) varies with similar ampli- 
tude (see Fig. 3). This indicates that either the value of k 
is indeed close to its lower limit, or the flux variations are 
governed by a different reason. In particular, one can also 
envisage a situation in which the magnetospheric radius of 
the neutron star during the low luminosity phase is close to 
the corotation radius (this implies k to be substantially larger 
than its minimum possible value) and the centrifugal barrier 
at the magnetospheric boundary reduces the mass transfer 
rate into the stellar magnetosphere (for discussion see Pema 
et al. 2006). The mass-transfer rate in the disk within this 
scenario varies insignificantly as the pulsar switches between 
the spin-up and spin-down phases. 

While the black body temperature during 2006-2010 in- 
creased from 023 keV to 0.51 keV, the source luminosity in- 
creased only by a factor of ~2. This implies that during the 
brightening phase the effective area of the black-body source 
has decreased by almost an order of magnitude. A high lumi- 
nosity and a relatively small size of this source suggest that 


5 



A. Camero-Arranz et al.: 4U 1626-67 as seen by Suzaku before and after the 2008 torque reversal 


it is located inside the magnetosphere in the vicinity of the 
neutron star. The accretion process in this region is fully con- 
trolled by the stellar magnetic field and plasma flows within 
the accretion channel. In this light, variations of the effec- 
tive area of the black-body source can be associated with a 
restructuring of the accretion channel (i.e. change of both 
its cross-section and the height of the shock at the base of 
the accretion column). Studies (see e.g. Bumard et al., 1983; 
Basko & Sunyaev, 1976; Kulkami & Romanova, 2008) sug- 
gest that the parameters of the channel depend on the ge- 
ometry and physical conditions in the accretion flow at the 
magnetospheric boundary, as well as on the mode by which 
the material enters the magnetic field of the neutron star. 
Our observations give no evidence of any significant change 
in the accretion flow geometry during the brightening phase 
(see above). It, therefore, appears more plausible to associate 
the restructuring of the accretion channel with change of the 
mode by which the accreting material enters the magneto- 
sphere. 

It is widely adopted that the accreting material can enter 
the magnetosphere due to the Rayleigh- Thy lor and Kelvin- 
Helmholtz instabilities or/and the magnetic reconnection. 
The Rayleigh- Taylor instability can be at work if the angular 
velocity of the accreting material at the boundary does not 
significantly differ from the angular velocity of the neutron 
star itself. The magnetic field configuration in the magne- 
topause in this case is close to poloidal. If this condition is not 
satisfied the Rayleigh-Taylor instability can be suppressed by 
the magnetic field shear in the magnetopause (Bumard et al., 
1983; Ikhsanov & Pustil’nik, 1996). This indicates that the 
relative velocity between the magnetic field of the neutron 
star and the accreting material at the boundary can be one 
cf the key parameters determining the structure of the ac- 
cretion channel. Even small fluctuations in the mass transfer 
rate may lead to significant variations of this parameter if the 
magnetospheric radius is close to the corotation radius (see 
e.g. Pema et al., 2006). Another possibility is an amplification 
cf the magnetic field at the inner radius of the disk governed 
by the dynamo action (see e.g. Rappaport et al., 2004). The 
disk in this case may switch into a sub-Keplerian state and 
decrease the angular velocity of the material at its inner ra- 
dius. However, the modeling of the evolution of the channel, 
and its influence on the properties of the black body source 
within these scenarios, is very complicated and beyond the 
scope of the present paper. 

£.2. Emission Lines 

The observed change of the equivalent width (EW) of the Ne Ly/3 
line at 1 .021 keV during the brightening phase of the pulsar from 
22.6 to 50.9 eV, as well as the increase of its intensity a factor of 
-8 reflect a transition of the accretion picture to a new phase. 
This may indicate either changes of plasma parameters at the 
inner radius of the accretion disk (Kii et al., 1986) or changes 
i;i the wind caused by interaction between the Keplerian disk 
and the stellar magnetic field at the corotation radius (Lovelace 
et a!., 1995). As shown by Pema et al. (2006), the mass outflow 
rate can be higher during the spin-down phase and lead to a pe- 
r.odic circulation of material between the inner and outer radii 
of the accretion disk. This feedback results in a redistribution of 
nass in the accretion disk and can lead to variations of the mass- 
transfer rate in the system. If this is the cause of the observed 
luminosity increase after the torque reversal in 2008, the time 
scale of radial motion of the material in the disk is r r ~ 2 years. 


Year 



Fig. 3. The X-ray flux history of 4U 1626-67 relative to the flux mea- 
sured by HE AO 7, in the same energy band, from previous works 
(Chakrabarty et al. (1997): circles; Orlandini et al. (1998): trian- 
gle; Krauss et al. (2007): stars; Camero-Arranz et al. (2011) unfilled 
squares) and the two Suzaku observations (blue diamonds) in the 0.7- 
60 keV band. 


This situation can be realized if the a parameter at the outer ra- 
dius of the disk is relatively small, 

a~0.02m 1/2 (— — T ) ’*(;£-) ^ < (2) 

\10 6 cms" 1 / \2yr/ V10 10 cm/ 

and, correspondingly, the density of material in the disk signifi- 
cantly exceeds its average value for X-ray pulsars (for discussion 
see e.g. King, 2012). Here c s is the sound speed in the disk at the 
distance R. The slow turbulization of the material in the disk can 
also explain the small value of the strength of the torque noise 
observed from 4U 1626-67 (Chakrabarty et al., 1997). 

In the present work, following Angeiini et al. (1995), we 
have decomposed the strong emission-line complex centered on 
1 keV found in the Suzaku X-ray spectrum of 4U 1626-67 (see 
Fig. 2), In both Suzaku observations the strongest line is the 
hydrogen-like Ne Lya line at -1 .025 keV, in agreement with the 
ASCA observations (Angeiini et al., 1995). The broad feature 
at ~1 .48 keV found in our observations is close to the expected 
energy of the H-like Ne recombination edge at 1.362keV. The 
measured 57 eV broadening of the line at ~1.48keV suggests 
that we may be also detecting the recombination continuum of 
Ne, as suggested by Angeiini et al. (1995) for the ASCA obser- 
vation. We point out that the Mg H a transition occurs at energy 
1.471 keV (which has not been resolved by Suzaku ), therefore 
this Mg line may be blended with the recombination contin- 
uum of Ne. Furthermore, with Suzaku an eighth line centered 
at ^1 .34keV was detected, which we tentatively identify as the 
Mg Hear (see Table 2). 

During the spin-down period, comparing the strength of the 
low-energy lines detected in the Mar 2006 Suzaku observation 
with previous measurements by ASCA , we see that the over- 
all decrease is of the order of ~2.5 (see Table 2). Before the 
2008 torque reversal, X-ray missions like BeppoSAX (Owens 
et al., 1997), Chandra (Schulz et al., 2001; Krauss et al., 2007) 
and XMM (Krauss et al., 2007) were able to resolve only up to 
3 or 4 lines of this complex. In general, the measurements of 
the intensity of those lines are in agreement with our results. 
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During the new spin-up period after the 2008 torque reversal, 
the flux of the lines increased in general by a factor of ~5 ex- 
cept the Ne Lya line at 1.021 keV that increased by a factor 
cf ~8. This result is not surprising taking into account such a 
change in luminosity. Unfortunately, during a previous spin-up 
period there was only one measurement of the ~1 keV complex 
by the Einstein observatory in 1979, and it was not resolved due 
to the lower data quality. The overall luminosity decrease asso- 
ciated with the 1990 torque reversal was a factor of -3 in the 
0.5-10keV range (Angelini et al., 1995, and references therein). 
Krauss et al . (2007) found a decrease of the equivalent widths 
of the NeHea 0.914keV and NeLya 1.021 keV lines during 
the spin-down period. The behavior of the EW for the OHeor 
0.568 keV and OLya 0.653 keV lines was more steady during 
the same period. With Suzaku we found a large increase of the 
EW of the Ne Lya 1 .021 keV line after the 2008 torque reversal. 
The EW for this line reached almost the same value measured 
by ASCA in 1993. For the other lines the increase is less pro- 
nounced. 

The strong emission-line complex at ~1 keV has been sim- 
ply explained as being due to a substantial Ne overabundance 
(Angelini et al., 1995). Angelini et al. (1995) raised the possibil- 
ity that 4U 1626-67 is embedded in a shock heated gas from the 
supernova remmant that created the pulsar, but as they point out, 
the absence of strong lines from Si , S, and Fe is inconsistent with 
this idea. Instead, Schulz et al. (2001) proposed that this cool ma- 
terial originated in the accretion disk or the mass donor, so that 
its composition is then a clue to the nature of the companion. 
Schulz et al. (2001) detected double-peaked emission lines and 
strong absorption edges of Ne and O in their Chandra HETGS 
spectrum, thus pointing to the companion in 4U 1626-67 being 
the chemically fractionated core of a C-0 or O-Ne-Mg crystal- 
lized white dwarf. Schulz et al. (2001) also concluded that the 
strong edges are due to absorption in cool, metal-rich material 
local to the source. Krauss et al. (2007) found that the double- 
peaked emission lines and the line ratios implied that they are 
formed in a high-density environment, suggesting that they arise 
somewhere in the accretion disk. Krauss et al. (2007) also noted 
that their value for the temperature of the line formation region 
fe 10 6 K) is characteristic of the highly ionized, optically thin 
outer layers of the accretion disk. We note that, in our work, sin- 
gle Gaussian line profiles provided an acceptable fit to the ob- 
served line shapes, and physical models like meka or vmeka did 
not yield good fits. 

Acknowledgements. A.CA. and M.H.F. acknowledge support from NASA 
grants NNX08AW06G and NNX1 1AE24G. D.M.M. and K.P. acknowledge sup- 
port from NASA grants NNX10AJ48G and NNX11AD41G N.R.I. acknowl- 
edges support from the Program N 21 of the Prezidium of RAS. 


Chakrabarty , D., Bildsten, L., Grunsfeld, J. M., et al. 1997, ApJ, 474, 414 

Ghosh, P. & Lamb, F. K. 1979, ApJ. 234, 296 

Giacconi, R, Murray, S., Gursky, H., et al. 1972, ApJ, 178, 281 

Ikhsanov, N. R. & Pustil’nik, L. A. 1996, A&A.312.338 

Iwakiri, W. B., Terada, Y., Mihara, T., et al. 2012, ApJ, 751 , 35 

Jain, C v Paul, B., & Dutta, A. 2010, MNRAS, 403, 920 

Kii, T., Hayakawa, S., Nagase, F., Ikcgami, T., & Kawai, N. 1986, PASJ,38, 751 

King, A. 2012, ArXiv e-prints 

Klochkov , D. , Doroshenko, V. , Santangelo, A. , et al. 2012, A&A, 542, L28 
Kommcrs, J. M., Chakrabarty, D. f & Lewin, W. H. G. 1998, ApJ, 497, L33+ 
Koyama, K. , T sunemi, H. , Dotani , T. , & et al. 2007, PAS J, 59, 23 
Krauss. M. I., Schulz, N. S., Chakrabarty, D., Juett, A. M., & Cottam, J. 2007, 
ApJ, 660, 605 

Kulkami, A. K. & Romanova. M. M. 2006, MNRAS, 386, 673 
Levine, A., Ma, C. P., McClintock, J., et al. 1988, ApJ, 327,732 
Lovelace, R. V. E., Romanova, M. M., & Bisnovatyi-Kogan, G. S. 1995, 
MNRAS, 275, 244 

McClintock, J. E., Brad:, H. V., Doxsey, R. E„ et al. 1977, Nature, 270, 320 
McClintock, J. E.,Li, F. K.,Canizares, C. R.,& Grindlay, J. E. 1980, ApJ, 235, 
L81 

Middlediteh. J., Mason, K. O., Nelson, J. E., & White, N. E. 1981, ApJ, 244, 
1001 

Mitsuda, K., Bautz, M.,Inoue, H.. & et al. 2007, PASJ, 59, 1 
Narayan, R. & Yi, I. 1995, ApJ, 452, 710 

Nowak, M. A., Hanke, M,, Trowbridge, S. N., et al. 2011 , ApJ, 728, 13 
Orlandini, M., Fiume, D. D., Frontcra, F, et al. 1998, ApJ, 500, L163+ 

Owens, A., Oosterbroek, T ., & Parmar, A. N. 1997, A&A, 324, L9 
Paczynski, B. & Sienkiewicz, R 1981, ApJ, 248, L27 
Pema, R. , Bozzo, E., & Stella, L. 2006, ApJ, 639, 363 
PoUschmidt, K., Suchy, S., Rivers, E., et al. 201 1, ArXiv e-prints 
Pravdo, S. H., White, N. E., Boldt, E. A., et al. 1979, ApJ, 231, 912 
Rappaport, S. A., Fregeau, J. M., & Spruit, H. 2004, ApJ, 606, 436 
Reynolds, A. P, Quaintrell, H., Still, M D., et al. 1997, MNRAS, 288, 43 
Schulz, N. S., Chakrabarty, D„ Marshall, H. L.,et al. 2001 , ApJ, 563, 941 
Shinoda, K., Kii, T., Mitsuda, K., et al. 1990, PASJ, 42, L27 
Shvaitsman, V. F. 1970, Radiophysics and Quantum Electronics, 13, 1428 
Suchy, S., Pottschmidi, K., Rothschild, R. E., et al. 2011, ApJ, 733, 15 
Takahashi, T., Abe, K., Endo, M.. & et al. 2007, PASJ, 59, 35 
Vaughan, B. A. & Kitamoto, S. 1997, ArXiv Astrophysics e-prints 
Wang, Y.-M. 1996, ApJ. 465, LI 11 
Wijers, R. A. M. J. & Pringle, J. E. 1999, MNRAS, 308, 207 
Wilson, R. B., Fishman, G. J., Finger, M. H., et al. 1993, in American Institute 
of Physics Conference Series, Vol. 280, American Institute of Physics 
Conference Series, ed. M. Friedlander, N. Gehreis, & D. J. Macomb, 291- 
302 

Yi, I. & Vishniac, E. T. 1999, ApJ. 516, L87 


References 

Angelini, L., White, N. E., Nagase, F, et al. 1995, ApJ, 449, L41 + 

Amaud, K. A. 1996, in Astronomical Society of the Pacific Conference 
Series, Vol. 101, Astronomical Data Analysis Software and Systems V, cd. 
G. H. Jacoby & J. Barnes, 17 
Basko, M. M. & Sunyaev, R A. 1976, MNRAS, 175, 395 
Becker, P. A. , Klochkov, D., Schonherr, G., et ai. 2012, ArXiv e-prints 
Bildsten, L., Chakrabarty, D., Chiu, J., et al. 1994, in American Institute of 
Physics Conference Series, Vol. 304, , 294-298 
Bozzo, E., Stella, L., Vietii, M„ & Ghosh, P. 2009, A&A, 493, 809 
Bumard.D. J., Arons, J.,& Lea, S. M. 1983, ApJ, 266, 175 
Caballero, I., Pottschmidt, K., Barragan, L., et al. 2009, in The Extreme Sky: 
Sampling the Universe above 10 keV 

Camero-Arranz, A., Finger, M. H., Ikhsanov, N. R, Wilson-Hodge, C. A., & 
Beklen, E. 2010, ApJ. 708, 1500 
Chakrabarty, D. 1998, ApJ, 492, 342 


7 



